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Abstract 


Land subsidence is one of the most important environmental problems that affects agriculture and 
urban infrastructure. In order to plan for mitigating the damages caused by land subsidence, it is 
necessary to identify high-risk and prone areas for this phenomenon specially in different plains 
of Iran. The main purpose of this research was zoning the land subsidence in Kashan aquifer using 
Dempster-Schafer, Weight of Evidence, Shannon Entropy, and Frequency Ratio. To conduct this 
research, 14 factors affecting the occurrence of subsidence including elevation, slope, slope 
direction, surface curvature, distance from stream, distance from fault, distance from road, 
lithology, land use, longitudinal and transverse curvature, water table changes, distance from 
mining and well density in aquifer surface were used. Then, determining 108 subsidence positions 
and 108 non-subsidence positions at the aquifer level, these points were randomly classified as 
30% and 70% as validation data and test data, respectively. To evaluate each of the models, the 
ROC curve was used. The results showed that among the used methods based on the ROC curve 
and the amount of AUC, the Frequency Ratio method (training: 0.84 and validation: 0.89) 
performed best and is introduced as the best method for predicting areas with subsidence 
sensitivity in studied area. Shannon Entropy, Dempster-Schafer and Weight of Evidence were 
also ranked next, respectively. The results of this study can be used by provincial managers to 
reduce the damage caused by land subsidence and to determine areas with subsidence sensitivity. 
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جغرافیا و مخاطرات محیطی. سال یازدهم. شمارة چهل و چهارم. زمستان ۱۶۰۱ صص 19-4۸ 
مقاله پژوهشی 
ویژه‌نامه (چالش جهانی فرونشست زمین: مدیریت بحران یا بحران مدیریت) 


نرگس سقازاده- دانش‌آموخته دکتری علوم و مهندسی آبخیزداری. دانشکده منابع طبیعی و علوم زمین, دانشگاه QUIS LEIS‏ ایران. 
o0‏ قاسمیه — دانشیار گروه مهندسی طبیعت» دانشکده منابع طبیعی 3 علوم CO)‏ دانشگاه کاشان» کاشان» ایران. 


ابراهیم امیدوار- دانشیار گروه مهندسی طبیعت. دانشکده منابع طبیعی و علوم زمین. دانشگاه کاشان کاشان ایران. 


پاسر مقصودی- دانشیار گروه فتوگرامتری و سنجش‌ازدور, دانشکده مهندسی نقشه‌برداری. دانشگاه خواجه‌نصیرالدین quo do‏ تهران ایران. 


تاریخ دریافت: ۱2۰۰/۱۱/۱۵ تاریخ بازنگری: ۱2۰۱/۲/۲۹ تاریخ تصویب: ۱۶۰۱/۳/۸ 


حکیده 


TU‏ زمین یکی از مهم‌ترین مشکلات محیط‌زیستی است که بر کشاورزی و زیرساخت‌های شهری 
تأثیر می گذارد. لذا به‌منظور برنامه‌ریزی برای کاهش خطرات ناشی از فرونشست. شناسایی مناطق پرخطر و 


مستعد ES‏ این پدیده در دشست‌های مختلف Ol pl‏ ضروری اسست. هدف T"‏ این تحفیق. 


زمین در آبخوان کاشان با استفاده از روش‌های دمپستر- شفر» وزن‌دهی 
شواهد. آنتروپی شانون و نسبت فراوانی است. برای انجام این تحقیق از ۱۶ عامل مؤثر بر وقوع فرونشست 
شامل ارتفا شیب جهت شیب انحنای سطح» فاصله از آبراهه» فاصله از گسل» فاصله از جاده. لیتولوژی» 
کاربری اراضی» انحنای طولی و عرضی. تغییرات سطح ایستابی, فاصله از معدن و تراکم چاه در سطح 
آبخوان استفاده شد. آنگاه پس از تعیین ۱۰۸ موقعیت فرونشست و نیز ۱۰۸ موقعیت عدم فرونشست در 
سطح آبخوان. این نقاط به‌صورت تصادفی و به نسبت ۳۰و ۷۰ درصد به‌ترتیب به‌عنوان داده‌های 


اعتبارسنجی و داده‌های آزمون طبقه‌بندی شدند و برای ارزیابی هرکدام از مدل‌ها از منحنی تشخیه 


Email: h.ghasemieh @kashanu.ac.ir ۰۹۱۳۲۳۱۳۷۷۹ نویسنده مسئول:‎ ۱ 

نحوه ارجاع به این مقاله: 

سقازاده, نرگس, قاسمیه, هدی, امیدوار, ابراهیم, مقصودی, یاسر. (۱۶۰۱). ارزیابی کارایی مدل-های دو متغیره در تعیین حساسیت 

پذیری فرونشست آبخوان دشت کاشان. جغرافیا و مخاطرات محیطی, (EV)‏ صص 1۹-۹۸ 
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Ve‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


عملکرد نسبی (ROC)‏ استفاده شد. نتایج نشان داد که از بین روش‌های مورداستفاده بر اساس منحنی 


می‌شود. روش‌های آنتروپی شانون» دمپستر- شفر و وزن‌دهی شواهد نیز به‌ترتیب در رتبه‌های بعدی قرار 


گرفتند. نتایج این تحقیق برای کاهش خسارات ناشی از وقوع فروند ت زمین و تعیین مناطق بالقوه 
فرونشست. LG‏ استفاده برای مدیران استانی خواهد بود. 


کلیدواژه‌ها: آب زیرزمینی» Cubs‏ کاشان» روش‌های دومتغیره. فرونشست. یاد گیری ماشینی 


۱- مقدمه 


ya‏ ت زمین ‏ که باعث تغییر شکا سطح زمین می‌شود پدیده خطرناکی است که می‌تواند خسارات بزرگی 


به زمین وارد کند (قربان‌زاده" و همکاران ۲۰۱۸). این پدیده ازجمله مهم‌ترین مخاطرات زمین‌شناسی است که 


در دهه‌های cue‏ گزارش‌هایی اعلام شده که تخلیه آبخوان بر کاهش AS‏ آب» فروند ت زمین و ترک 


خوردن زمین SE‏ گذاشته است (ژو و لی“ ۲۰۱۱). فرونشست زمین یک پدیده جهانی است که از فرایندهای طبیعی 
مانند تراکم خودکار و فعالیت‌های زمین‌ساختی و همچنین فعالیت‌های انسانی که منجر به خروج و بهره‌برداری 
بی‌رویه از آب» CH‏ و سایر سیالات از درون زمین می‌شود ناشی می‌شود و می‌تواند آسیب جدی به زیرساخت‌های 
اقتصادی بزند (گالوی* و همکاران. ۲۰۱۳ SU‏ و همکاران ۲۰۲۲). 

به نظر می‌رسد که در بین عوامل مؤثر در ایجاد فرونشست. برداشت بیش‌ازحد مجاز از منابع آب زیرزمینی عامل 
اصلی ایجاد فرونشست در بیشتر دشت‌های ایران باشد؛ به‌گونه‌ای که در بیست و یکمین کارگروه ملی فرونشست 


زمین و حفاظت از منابع آب زیرزمینی, بیشتر دشت‌های کشور در معرض خطر فرونشست بر اثر برداشت بی‌رویه از 


ples‏ آب زیرزمینی تشخیص داده شدند (رنجبر و جعفری» ۸ طی پژوهش‌های صورت گرفته توسط محققین. 


بهره‌برداری و استخراج آب‌های زیرزمینی در چندین کشور در سراسر زمین شامل مکزیک. Ol pl Wh!‏ ویتنام و 


1 Land Subsidence 
2 Ghorbanzadeh 

3 Ashraf & Cawood 
4 Zhou & Li 

5 Galloway 

6 Liang 


سال eh‏ ارزیابی کارایی مدل‌های دو متغیره در Oe‏ حساسیت پذیری — VA]‏ 


چین» نقش مهمی در ایجاد فرونشست زمین oa S Lal‏ است e)‏ ' و همکاران» ۲۰۲۰؛ ژو" و همکاران ۲۰۲۰؛ تای " و 
همکاران» ۲۰۲۱؛ oy gi‏ و همکاران» ۲۰۲۲. اکبری» ۱۶۰۰؛ قرهچلو و همکاران ۱۶۰۰). 

در سال‌های tl‏ روش‌های یادگیری ماشینی * و روش‌های دومتغیره" همراه با استفاده از سامانه اطلاعات 
جغرافیایی در مقایسه با روش‌های سنتی در زمینه تعیین مناطق حساس به فرونشست که poles‏ هستند. دارای کارایی 
بالایی بوده‌اند و باعث افزایش دقت و سرعت در زمینه ارزیابی و پتانسیل فرونشست شده‌اند. فاکتورهای مختلفی از 
طیف وسیعی از پارامترهای محیطی از قبیل فاکتورهای هیدرولوژیکی. gl‏ اقلیمی, فاکتورهای توپوگرافیک» 
clay SU‏ اکولوژیک فاکتورهای ژئولوژیکی و فاکتورهای انسانی بر میزان حساسیت فرونشست S SU‏ هستند. 
در این dhe‏ روش‌های یادگیری ماشینی (مدل‌های جنگل تصادفی " حداکثر oy I‏ شبکه عصبی مصنوعی* و 
ماشین بردار پشتیبان ") و دومتغیره (مدل‌های نسبت فراوانی ٩‏ دمیستر- شفر c‏ وزن‌دهی شواهد " و رگرسیون 
خطی *' و لجستیک*؟) وجود دارد که به محققان اجازه می‌دهد نقشه حساسیت فرونشست را تهیه کنند و با تعیین 
روابط آماری بین عوامل محیطی و توزیع مکانی فرونشست‌هاء احتمال مکانی وقوع فرونشست را ارزیابی و پیش‌بینی 
aus‏ در سال‌های اخیر از تکنیک‌های رایج ازجمله یادگیری ماشین و روش‌های آماری و دو متغیره برای تخمین 
میزان فرونشست و بررسی علل این پدیده استفاده شده است. به‌عنوان‌مثال پرادهان و ابوخاریما"" OVE)‏ در تحقیقی, 


میزان فرونشست زمین را در دره کینتاپراک مالزی با استفاده از GIS‏ و روش‌های سنجش ازدور برآورد نمودند. در این 


تحقیق» ۱۶ شاخحص به‌عنوان عوامل $$ بر میزان 33 ت در نظر گرفته شد. آن‌ها به‌منظور تهیه نقشه 
حساسیت‌پذیری فرونشست از دو روش تابع باور شهودی" و نسبت فراوانی استفاده کردند و ارزیابی مدل‌ها در این 
پژوهش با استفاده از سطح زیر منحنی (AUC) “ROC‏ اعتبارسنجی شد. یافته‌های این پژوهش نشان داد که مدل 


1 Yu 

2 Zhu 

3 Ty 

4 Nguyen 

5 Machine Learning 

6 Bivariate 

7 RF: Random Forest 

8 MF: Maximum Entropy 

9 ANN: Artificial Neural Network 

10 SVM: Support Vector Machine 

11 FR: Frequency Ratio 

12 Dempster-Shafer 

13 WOF: Weight Of Evidence 

14 Linear regression 

15 Logistic Regression 

16 Pradhan & Abokharima 

17 EBF: Evidential Belief Function 
18 ROC: Relative Operating Curve 
19 AUC: Area Under Curve 
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نسبت فراوانی قادر به ایجاد نرخ پیش‌بینی ۷۱/۱۲ درصد است. درحالی‌که مدل تابع باور شهودی پیش‌بینی بهتری را 
با نرخ ۷۳/۱۳ درصد نشان داد. تین‌بویی " و همکاران (۲۰۱۸) با استفاده از چهار مدل رگرسیون لجستیک بیزین " 
ماشین بردار پشتیبان درخت مدل لجستیک y‏ درخت تصمیم متناوب y‏ ۸ عامل مور حساسیت فرونشست زمین 
در کره جنوبی را ارزیابی کردند و به این نتیجه رسیدند که مدل رگرسیون لجستیک بیزین از دقت بالاتری برخوردار 
بوده است. همچنین برداشت بی‌رویه آب زیرزمینی را عامل فرونشست در منطقه موردمطالعه بیان کردند. رحمتی* و 
همکاران (۲۰۱۹۵) به بررسی خطر فرونشست در دشت همدان با استفاده از چهار مدل درخحت تصمیم قانون محور 
درخت رگرسیون تقویت شده" رگرسیون درختی و طبقه‌بندی" و جنگل تصادفی پرداختند. آن‌ها برای انجام 
مدل‌سازی و ag‏ نقشه das‏ فرونشست از شش عامل استفاده کردند. درنهایت بیان داشتند که روش جنگل تصادفی 
بهترین نتیجه را داشته و کاهش سطح آب زیرزمینی, لیتولوژی و فاصله از آبراهه‌ها مهم‌ترین عوامل موثر بر فرونشست 
iles‏ 

caa‏ دشک ed dl pts‏ بارس مال ۱۳۷۶ uo deg‏ شکاف‌هایی در duis 3l alu es‏ متاح 
گردید که باعث تخریب واحدهای مسکونی و زمین‌های کشاورزی گردید (قاضی‌فرد و همکاران» ۱۳۹۰). علت این 
عوارض در آن زمان. استخراج بی‌رویه آب زیرزمینی عنوان شد. اما به دلیل گسترش شکاف‌ها و پیشروی آن‌ها در 
بسیاری از نقاط دشست. بررسی جامع این پدیده به‌منظور شناخت JAS‏ این عوارض,. تعیین محدوده تأثیرپذیر 
پیش‌بینی نحوه گسترش و در صورت امکان» جلوگیری از توسعه و تشدید آن‌هاء پيشنهاد گردید (مصلحی بهارانچی 
DUCI uae à‏ ور ca‏ خی لاش هی شوه با قاد از sita gs cd Unda‏ سای ات pts‏ 


در آبخوان دشت کاشان ag‏ شود تا مناطق das s‏ شناسایی شوند. در نظر گرفتن ۱۶ عامل مؤثر و به‌کارگیری آن‌ها 


در مدل‌های دو متغیره برای تهیه نقشه‌های حساسیت‌پذیری فرو در آبخوان دشت کاشان از نوآوری‌های 


تحقیق حاضر محسوب می‌شود. 


] Tien Bui 

2 BLR: Bayesian Logistic Regression 

3 LMT: Logistic Model Tree 

4 ADT: Alternate Decision Tree 

5 Rahmati 

6 RBDT: Rule-Based Decision Tree 

7 BRT: Boosted Regression Trees 

8 CART: Classification And Regression Tree 


Y‏ مواد و روش 

۱-۲- منطقه موردمطالعه 

حوزه آبخیز دشت DUIS‏ در جنوب حوزه آبخیز بزرگ دریاچه نمک بین طول جغرافیایی ۵۰ درجه و ۵۶ دقيقه 
تا ۵۲ درجه و " دقيقه شرقی و عرض جغرافیایی ۲۳ درجه و ۳۱ دقیقه تا ۳۶ درجه و ۳۰ aiia‏ شمالی واقع شده 
است. این حوضه دارای وسعت 1۹۵۰ کیلومترمربع است. دشت کاشان در دامنه کوه‌های کرکس و حاشیه US‏ 
مرکزی olal‏ در حدود ۲۶۰ کیلومتری جنوب تهران و بین طول‌های جغرافیایی ۵۱ و ۵۲ درجه شرقی و عرض‌های 
۰ و ۳۶/۳۰ درجه شمالی درجه قرار دارد (شکل ۱). آبخوان دشت LES‏ دارای وسعت ۲۲۱۷۳۱۳/۶ هکتار است 
که البته بخش‌هایی از آن به دلیل شور بودن آب‌های زیرزمینی» قابل‌استفاده نیست. دشت کاشان, به‌صورت یک دره 
باریک دارای امتداد شمال غربی- جنوب شرقی است. حداقل و حداکثر ارتفاع دشت OUS‏ از سطح دریاء به‌ترتیب 
۰و ۱۲۰۰ متر در قسمت حاشیه شرقی و غربی دشت کاشان است (میرزاوند و OL [Kar‏ ۱۳۹۷). دامنه تغییرات 
ارتفاعی منطقه مطالعاتی در این پژوهش نیز بین VWI EASY.‏ متر است و از شمال و شمال غربی منطقه مطالعاتی به 
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شکل ۱- موقعیت دشت کاشان در استان اصفهان و ایران 


Ví‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


۲-۲- روش‌شناسی GAPS‏ 
۱-۲-۲- تهیه نقشه پراکنش نقاط وقوع و عدم وقوع فرونشست 
ابتدا Modum az T po‏ تقاط وقوع 3 عدم ES‏ فرونشست. بازدیدهای میدانی از سطح دشت به‌ویژه در 


مناطقی که چاه‌های پیزومتری و بهره‌برداری وجود داشت. انجام گرفت و درنهایت تعداد ۸ مورد 99 23 


مشخص شد که با استفاده از GPS‏ موقعیت جغرافیایی این نقاط برداشت و ثبت شد. به‌منظور ارزیابی مدل‌های دو 
متغیره علاوه بر این ۱۰۸ مورد تعداد ۱۰۸ نقطه فاقد فرونشست نیز طی عملیات صحرایی ثبت شد (شکل (Y‏ ذکر 
این نکته ضروری است که به‌منظور انتخاب اولیه این نقاط ابتدا از نقشه میزان فرونشست دشت که در مطالعه‌ای 
موازی به روش تداخل‌سنجی راداری توسط نویسندگان تهیه شده استفاده گردید و مناطق فاقد فرونشست شناسایی 
شدند. سپس طی عملیات صحرایی نقاطی که فاقد هرگونه آثار فرونشست در اطراف خود بودند به‌عنوان نقاط عدم 
iy‏ 

مخروط افکنه‌های Gl Dl‏ شجاع‌آباد و احمدآباد (فاصله بین جاده قدیم بادرود و اتوبان اصلی) است که فرونشست در 


ت ثبت شدند. لازم به ذکر است که بخش جنوب غربی آبخوان دشت کاشان همان منطقه تپه‌ماهوری و 
آن‌ها مشاهده نشد. 
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شکل ۲- نقشه پراکنش نقاط وفوع و عدم وقوع فرونشست 


سال eh‏ ارزیابی کارایی مدل‌های دو متغیره در Oe‏ حساسیت پذیری یک ۷۵ 


ang -۲-۲-۲‏ نقشه عوامل مؤثر بر وقوع فرونشست 

مطابق با سوابق پژوهشی (مختاری و همکاران. ۱۳۹۸؛ محمدی" و همکاران» ۲۰۱۹ قربان‌زاده و همکاران ۲۰۱۸؛ 
عبداللهی و همکاران» ۲۰۱۹؛ رحمتی " و همکاران. (Y 19D‏ و بازدیدهای میدانی به‌منظور انجام مدل‌سازی ميزان 
فرونشسست در سطح آبخوان دشت کاشان» تعداد VE‏ عامل مؤثر بر وقوع فرونشست در چهار گروه کلی عوامل 
توپوگرافیک (شیب. جهت. ارتفاع از سطح دریا انحنای طولی دامنه " انحنای عرضی aal‏ و انحنای ترکیبی دامنه )» 
عوامل زمین‌شناسی (سنگ‌شناسی و فاصله از گسل)» عوامل انسانی (کاربری اراضی فاصله از جاده تراکم چاه‌های 
بهره‌برداری و فاصله از معادن) و همچنین عوامل هیدرولوژیک (تغییرات سطح ایستابی و فاصله از آبراهه) در نظر 


گرفته شدند. 


- عوامل توپوگرافیک 

به‌منظور تهیه نقشه عوامل توپوگرافیک از نقشه‌های مدل رقومی ارتفاع آبخوان دشت کاشان با قدرت تفکیک 
مکانی ۱۳/۵ متر استفاده شد. این نقشه‌ها مستخرج از داده‌های سنجش ‌ازدور ماهواره MALOS‏ بوده و از پایگاه 
https://vertex.daacc.asf.alaska.edu‏ دانلود شد. یکی از عوامل مؤثر بر وقوع فرونشست. عامل ارتفاع از سطح Ly‏ 
است. تغییرات ارتفاع موجب تغییرات بارندگی می‌شود و درنتیجه بر میزان هوازدگی و میزان نفوذ آب تأثیرگذار است 
(تیموری و اسدی نیلوان» YA‏ رحمتی و همکاران ۲۰۱۹2). این نقشه با استفاده از دستور Reclassify‏ روی نقشه 
مدل رقومی ارتفاع در پنج طبقه ارتفاعی به دست آمد. عامل شیب نیز یکی از متغیرهای پستس و بلندی خیلی مهم در 
وقوع فرونشست زمین است (رحمتی و همکاران ۲۰۱۹۵). جهت شیب نیز بر میزان دریافت و نفوذ بارند گی» تبخیر 
و تعرق و میزان رطوبت اثرگذار است و بر همین اساس, نوع پوشش گیاهی و توزیع آن را تحت SE‏ قرار می‌دهد 
(شفاپور تهرانی " و همکاران» ۲۰۱۷). لایه‌های شیب و جهت شیب با استفاده از مدل رقومی ارتفاعی مورداشاره و به 
ترتیب از Slope els‏ و Aspect‏ در محیط نرم‌افزار 10.4 ARCGIS‏ تهیه شدند. نقشه‌های انحنای عرضیی انحنای 
Job‏ و انحنای عمومی از مدل رقومی ارتفاعی دشت و با استفاده از دستورهای مربوطه در محیط ang Arc GIS‏ 
شدند. در این مطالعه برای طبقه‌بندی نقشه‌های این عوامل از روش شکست طبیعی ^ استفاده شد. این شاخحص‌ها در 


1 ۷ 

2 Rahmati 

3 Profile Curvature 

4 Plan Curvature 

5 Curvature 

6 Advanced Land Observation Satellite 
7 Shafapour Tehrany 

8 Natural break 


vi‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


باعث تجمع آب و در شیب‌های محدب باعث افزايش رواناب می‌شوند (نقیبی و همکاران؛ ۲۰۱۶؛ تیموری و اسدی 


(ال- ابادی و همکاران؛ ۲۰۱7؛ عبداللهی و همکاران ۲۰۱۹). 


- عوامل زمین‌شناسی 


نوع سنگ و سازندهای موجود در یک منطقه» عامل بسیار مهمی در وقوع فرونشست محسوب می‌شوند و 
واحدهای سنگی مختلف دارای درجات حساسیت گوناگونی در وقوع فرونشست زمین هستند (رحمتی و همکاران؛ 
“OL SES ۹ a‏ و همکاران. ۲۰۱۶؛ زمانی‌راد؛ و همکاران. ۲۰۲۰). در پژوهش حاضر aU‏ سنگ‌شناسی منطقه 
موردمطالعه از نقشه‌های زمین‌شناسی شیت‌های OLE‏ آران و کهک با مقیاس ۱:۱۰۰۰۰۰ که توسط سازمان 
زمین‌شناسی کشور تهیه شده استخراج گردید. وجود گسل نیز تأثیر زیادی در فرونشست زمین دارد (قربان‌زاده و 
همکاران» ۲۰۱۸؛ رحمتی و همکاران. ۲۰۱۹0 مختاری و همکاران» ۱۳۹۸). بدین‌منظور از عامل فاصله از گسل استفاده 
شد. با توجه به این که برخی از گسل‌ها در فاصله نزدیکی خارج از منطقه مورد قرار دارند و می‌توانند تأثیر خود را در 
نواحی حاشیه‌ای بگذارند؛ لذا در این پژوهش. گسل‌های اطراف منطقه نیز موردبررسی و توجه قرار گرفتند. همچنین 
در روش‌های دو متغیره مورداستفاده در این تحقیق, محاسبات بر اساس طبقه‌بندی عوامل موثر انجام گردید. ازاین‌ری 
این متغیر نیز به‌صورت طبقات مختلف فاصله از گسل ارائه گردید و درنهایت ada‏ فاصله از گسل با توجه به محل 
slo LS‏ اصلی و فرعی داخل و اطراف دشت و با استفاده از دستور Multi ring buffer‏ با فواصل ۵۰۰ متری در 


نرم‌افزار 10.4 Arc GIS‏ تهیه شد. 

- عوامل انسانی 

ازجمله عوامل انسانی موثر در وقوع فرونشسست. کاربری اراضی است. این عامل روی رواناب سطحی و میزان 
نفوذ آب اثرگذار است و از طریق مقدار نفوذ بر وقوع فرونشست تأثیر می‌گذارد. درنتیجه» پوشش ALS‏ نقش مهمی 
در پایداری دشت دارد (رحمتی و همکاران ۲۰۱۹ در cul‏ تحقیق برای عامل کاربری اراضی از ALAS‏ تهیه شده 
برای مطالعات جامع حوضه کاشان که توسط اداره کل منابع طبیعی و آبخیزداری استان اصفهان در سال ۱۳۹۳ تهیه 
شده oy‏ استفاده گردید. همچنین در این مطالعه. قبل از به‌کارگیری» نقشه مذکور به فرمت kml‏ تبدیل و طی 
هم‌پوشانی با تصویر متناظر گوگل ارث منطقه. به‌صورت بصری از دقت و صحت نسبی آن اطمینان حاصل گردیده 


است. 


1 Naghibi et al. 
2 Al-Abadi 

3 Kaitantzian 
4 Zamanirad 


سال eh‏ ارزیابی کارایی مدل‌های دو متغیره در Oe‏ حساسیت پذیری VY Meus‏ 


فاصله از جاده نیز یکی دیگر از شاحص‌های انسانی است که اثرات آن روی تغییرات سطح زمین می‌تواند منجر 
به فرونشست شود. در این مطالعه. نقشه شبکه جاده‌های داخل و خارج از دشت نیز از تصاویر Google Earth‏ تهیه و 
همانند نقشه‌های فاصله از آبراهه‌ها و گسل‌ها با استفاده از ابزار Multi ring buffer‏ به نقشه فاصله از جاده با فواصل 
۰ متری des‏ شدند. 

تعداد بیشتر چاه نیز باعث افزایش تراکم چاه و درنتیجه برداشت بیشتر از آبخوان» به کاهش سطح آب زیرزمینی 
آبخوان منجر خواهد شد (قربان‌زاده و همکاران» ۲۰۱۸). در این راستاء aU‏ تراکم چاه‌ها بر اساس چاه‌های بهره‌برداری 
موجود در منطقه و با استفاده از تابع Point Density‏ بر حسب تعداد در کبلومترمربع در نرم‌افزار 104 ag Arc GIS‏ 

دیگر عامل انسانی» احداث معدن است که از طریق لرزش‌های ناشی از انفجارات معدن‌کاوی» می‌تواند روی 
فرونشست دشت تأثیر گذار باشد و منجر به از دست رفتن مواد مغذی خاک و کاهش حاصلخیزی خاک شود (سان! 
و همکاران. ۲۰۲۳). در این تحقیق, این عامل به‌صورت a‏ فاصله از معدن در نظر گرفته شده و نقشه آن با استفاده از 
موقعیت نقاط معادن داحل و خارج از محدوده دشت و با استفاده از دستور Multi ring buffer‏ با فواصل مساوی 


۰ متری در نرم‌افزار 10.4 Arc GIS‏ تهیه شد. 


- عوامل هیدرولوژیک 

تغییرات سطح ایستابی آب زیرزمینی یکی از مهم‌ترین عوامل مؤثر بر فرونشست است. در هنگام افت سطح آب 
e‏ زهکشی col‏ سبب کاهش فشار هیدرواستاتیکی در all‏ رس و تحکیم و فشردگی این لایه‌ها و نشست 
زمین می‌شود (کیانی» (YAT‏ به‌منظور تهیه نقشه تغییرات سطح ایستابی, ابتدا آمار سطح ایستابی چاه‌های پیزومتری 
واقع در منطقه برای یک دوره ۱۲ ساله (۱۳۸۶-۹۲) از مدیریت منابع آب کاشان اخذ و سپس در هر نقطه چاه با کسر 
رقوم سطح آب در سال ٩۱‏ از رقوم سطح آب در سال ۸۶ میزان تغییرات سطح ایستابی در هر چاه به دست آمد. 
درنهایت با استفاده از روش میان‌یابی وزن‌دهی معکوس فاصله TDW)‏ نقشه تغییرات سطح ایستابی تهیه شد. 


موقعیت آبراهه‌ها نیز به‌عنوان یک عامل هیدرولوژیک» نقش بسیار مهمی در فرو زمین ایفا می LS‏ 
(قربان‌زاده و همکاران ۲۰۱۸؛ رحمتی و همکاران, b‏ ٩۲۰۱؛‏ مختاری و همکاران ۱۳۹۸). آبراهه‌ها می‌توانند نقش 
تغذیه سفره‌های آب زیرزمینی را در مواقع سیلابی ایفا نماینده ازاین‌رو مناطق نزدیک به آبراهه می‌توانند از این مزیت 
برخوردار شوند. در این راستا از عامل فاصله از آبراهه استفاده شد. برای تهیه نقشه فاصله از آبراهه ابتدا لازم بود نقشه 


شبکه زهکشی منطقه موردمطالعه تهیه شود. بدین‌منظور از ابزار Arc Hydro‏ استفاده و نقشه مذکور تهیه شد. در ادامه, 


1 Sun 
2 IDW: Inverse Distance Weighting 


VA‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


ایه‌های فاصله از آبراهه بر اساس نقشه آبراهه‌ها و با استفاده از دستور Multi ring buffer‏ با فواصل ۲۰۰ متری در 
نرم‌افزار ag Arc GIS‏ شد. لازم به ذکر است که با مرور ale‏ و سوابق تحقیقاتی مشخص LoS‏ که با توجه به 
وسعت و تراکم عوارض خطی در هر منطقه فواصل در مقالات مختلف متفاوت در نظر گرفته شده است. لذا در این 
پژوهش چندین و چند حالت برای فواصل در نظر گرفته شد و فواصل بافرهای عوارض خطی بر اساس یک فاصله 
بهینه انتخاب شد تا طبقات با فرونشست صفر حداقل ممکن شود. درنهایت این فواصل بر اساس واقعیت و شرایط 
منطقه انتخعاب شدند. 

۳-۲- تعیین حساسیت‌پذیری فرونشست با استفاده از مدل‌های دو متغیره 

در این مرحله از پژوهش, از مجموع ۱۰۸ موقعیت فرونشست ثبت شده به‌صورت تصادفی VA‏ نقطه فرونشست 


(۷۰ درصد) به‌عنوان داده‌های آموزش برای فرایند مدل‌سازی و 059( loa‏ انتخاب شدند. همچنین تعداد YY‏ 


نقطه فرونشست (۳۰ درصد) باقی‌مانده نیز برای اعتبارسنجی (آزمون) مدل‌ها در نظر گرفته شدند (شکل ۲). سپس 
به‌منظور وزن‌دهی به هر یک از طبقات عوامل مؤثر ۱۶ گانه. از مدل‌های دو متغیره دمپستر- شفر (تابع باور شهودی) 
آنتروپی شانون نسبت فراوانی و وزن شواهد استفاده شد که هر یک از آن‌ها در ذیل تشریح می‌شود: 

۱-۳-۲- مدل دمپستر- شفر (تابع باور شهودی) 

تئوری دمپستر- شفر به‌عنوان تعمیمی از تئوری بیزین است که مبتنی بر احتمالات است. این تئوری در سال ۱۹٩۷‏ 
توسط دمپستر به کار برده شد (دمپستر NAW‏ صابرچناری و همکاران» (YAT‏ و در ادامه» ساختار ریاضی آن توسط 
شفر در سال ۱۹۷ تشریح گردید (شفر V ٩‏ $10 فرض شود چندین عامل از داده‌های مکانی موجود در منطقه 


حساس به فرو است» هر لایه‌ای از داده‌های مکانی به‌صورت )1=1,2,......1( Ei‏ برای تابع هدف Tp‏ در نظر 


گرفته می‌شود. به‌این‌ترتیب Ej‏ به‌دست می‌آید که در آن زویژگی کلاس یا طبقه‌ای از Ei‏ دارای توزیع فراوانی از 
موقعیت‌ها L‏ ویژگی‌های مت و منفی در بروز فرونشست‌هاست. c»‏ مثبت با علامت XCTp)g;‏ نشان داده می‌شود 
که از رابطه ۱ به‌دست می‌آید (پارک " QV‏ 


= N(L) 
و(‎ X N(Ejj)-N(LNEjj) Q) 
NANA 


نرخ یا نمره منفی نیز با علامت ATP) E‏ نیز از رابطه ۲ به‌دست می‌آید: 


] Dempster 
2 Shafer 
3 Park 


سال يازدهم ارزیابی کارایی مدل‌های دو متغیره در Oe‏ حساسیت پذیری Meus‏ ۷۹ 


N(L)-N(LNE;;) 


7 = N(L) 
AT), ^. N()-NQ)-N(Eij)*N(LnE;;) (۳ 
N(A)-N(L) 
«iat iere 225 À عواملی است که بر روی‎ Tp «c 9 عوامل موثر در ایجاد‎ Tp که در آن‎ 


NÒL N Ey)‏ مساحت فرونشست در هر طبقه» N(L)‏ مساحت کل فرونشست‌های منطقه. (Ej)‏ مساحت فاقد 
فرونشست در هر طبقه و N(A)‏ مساحت فاقد فرونشست منطقه مطالعاتی است. به‌طورکلی مقدار نرخ‌های به‌دست 


آمده از رابطه‌های ۱ و Y‏ در محدوده صفر تا بی‌نهایت متغیر است. لدا با استفاده از رابطه‌های ۳ € و ۵ استانداردسازی 


صورت می‌گیرد: 
,رو( ATP‏ 
uj‏ = 
m(Tp)z., = DACP )E,, (v)‏ 
A(p)g,;‏ — 
a‏ 
)£( و۵( 207 2۱ mm‏ روم 1701 


در این رابطه‌هاء m(Tp)gy‏ میزان اعتقاد به وقوع فرونشست: m(Tp)g,‏ میزان اعتقاد به عدم وقوع فرونشست و 
m(9)‏ میزان نرخ نهایی هستند (پارک» ۲۰۱۰). 

۲-۳-۲- روش آنتروپی شانون 

روش آنترویی شانون در سال ۱۹۶۸ bug‏ شانون ارائه شد. این روش به‌عنوان یکی از معروف‌ترین روش‌های 
تصمیم گیری چند معیاره برای محاسبه وزن‌های معیارها و شاخص‌ها شناخته شده است. در این روش معمولاً به هر 
یک از شاخص‌ها یک وزن نسبت داده می‌شود به‌طوری‌که مجموع وزن شاخص‌ها برابر ۱ باشد. در این روش هر 
چه پراکندگی در مقادیر یک شاخحص بیشتر Lab‏ آن شاخحص از اهمیت بیشتری برخوردار است. شانون نشان داد که 
وقایع با احتمال وقوع ob‏ اطلاعات کمتری در اختیار می‌گذارند و برعکس هرچقدر احتمال وقوع یک رخداد کمتر 
باشب اطلاعات حاصل از Of‏ بیشتر است. با به‌دست آوردن اطلاعات جدید در واقع عدم قطعیت‌ها کاهش یافته و 
ارزش اطلاعات جدید برابر با مقداری است که از عدم قطعیت کاسته شده است. درنتیجه. عدم قطعیت y‏ اطلاعات 
پارامترهایی وابسته به هم هستند. در این پژوهش به‌منظور تعبین میزان حساسیت فرونشست و ترسیم نقشه با استفاده 
از مدل آنتروپی شانون از رابطه‌های زیر استفاده شد (یورقاسمی و همکاران ۲۰۱۲): 


(Pj) = Pi/ Xj, Pj C) 


^ جغرافیا و مخاطرات محیطی شمارة چهارم 


Himax = log» رگ‎ (A) 
| = e p -(01. j=12..n (4) 
Hjmax 


که در Pi iod‏ تراکم پا نسبت فراوانی (FR)‏ فرونشست در هر طبقه است که در رابطه ۱۰ نحوه محاسبه آن توضیح 


داده می‌شود. Pij‏ احتمال وجود فروند ت در هر عامل و کلاس‌های مرتبط با آن Hj‏ و Him‏ مقادیر آنتروپی و 
ماکزیمم آنتروپی؛ 1 ضریب اطلاعات؛ رک تعداد طبقه‌ها و Wj‏ وزن نهایی هر عامل است. بعد از تعیین وزن نهایی هر 
عامل و ضرب آن در کلاس‌های عامل مذکور یعنی مقدارهای Pi‏ مربوط به هر عامل و کلاس‌های آن بر اساس تقسیم 
تعداد فرونشست بر تعداد فرونشست بر تعداد پیکس‌های کلاس‌های عوامل» نقشه‌های وزنی با هم جمع شده و نقشه 
نهایی پتانسیل فرونشست تهیه شد (ذبیحی و همکاران, ۱۳۹۶). 

۳-۳-۲- روش نسبت فراوانی (FR)‏ 

نسبت فراوانی یک روش آماری دو متغیره است که برای محاسبه رابطه احتمالی بین متغیرهای مستقل و وابسته 
در نقشه‌های با کلاس‌های متعدد استفاده می‌شود و می‌تواند روابط بین فرونشست‌ها y‏ عوامل مختلف مؤثر در منطقه 
موردمطالعه را آشکار نماید (ازدمیر و آلتورال؛ ۲۰۱۳). در این پژوهش برای محاسبه نسبت فراوانی هر یک از عوامل 
مؤثر بر وقوع فرونشست. درصد پیکسل‌های دارای فرونشست و فاقد فرونشست منطقه موردمطالعه با استفاده از 


سامانه اطلاعات جغرافیایی تعیین گردید. سپس نسبت فراوانی برای هر یک از عوامل و کلاس مرتبط با آن محاسبه 
TS) da $‏ ۲۰۰۷؛ محبی تفرشی و همکاران, 09۹ محاسبه نسبت فراوانی برای هر یک از کلاس‌های عوامل مور 


A 
m 28 0۱0 


مطالعاتی» IC‏ مساحت هر کلاس مربوط به هر عامل و D‏ مساحت کل منطقه موردمطالعه است. به‌عبارت‌دیگر» نسبت 
بین پهنه‌های دارای فرونشست به مساحت کل منطقه بیان‌گر نسبت فراوانی است. ازاین‌رو Fr<l‏ نشان‌دهنده 
همبستگی ضعیف بین فرونشست‌های رخ‌داده و عوامل موثر بر وقوع آن و 1< Fr‏ نشان‌دهنده همبستگی زیاد بین 
فرونشست‌های رخ‌داده و عوامل موثر بر وقوع آن است. Fr‏ نزدیک به یک نشان می‌دهد که همبستگی معنی‌داری 


1 Ozdemir 6 ۷ 
2 Lee 


سال يازدهم ارزیابی کارایی مدل‌های دو متغیره در Oe‏ حساسیت پذیری — ^Y‏ 


بين 3 233 ت‌های رخداده و عوامل موثر بر وقوع Ol‏ وجود ندارد. به‌طورکلی» همبستگی مثبت نشان‌دهنده احتمال 
TC‏ متا بالاتر و ^ همبستگر منفی نشان‌دهنده احتمال A 3j‏ تست پایین‌تر است Yan «SU ple)‏ پرادهان و لی 4 
Qi‏ 


-£-Y-Y‏ مدل وزنی شاهد (وزن شواهد) 


این روش برای مدل‌سازی در یک منطقه برای هر عامل در صورت وجود يا عدم وجود فرونشست در محل» 
وزنی را محاسبه می‌کند. تعیین وزن برای هر عامل از طریق ترکیب هرکدام از لایه‌های عاملی با نقشه پراکنش 


انجام می‌شود. در این مدل» مقادیر وزنی معمولاً دامنه‌ای بین اعداد مثبت و منفی دارند که بیان‌گر نقش 


ua‏ (تویید " و همکاران. ۲۰۰۷؛ علیخانی و همکاران, ۰ وزن‌ها در این روش بر اساس وجود يا عدم وجود 


ü P 
w* = Ine) (11) 


P(— 
wW- =n (VY) 


C=W*-W" aY) 


که در آن P‏ احتمال, ۴ حضور یک فاکتور: *17عدم حضور یک فاکتور S‏ حضور فرونشست و S*‏ عدم حضور 
فرونشست است. Wh‏ و ۷۷ به ترتیب بیان‌گر وزن شواهد وقتی یک فاکتور (مثلاً زمین‌شناسی)» حضور و عدم حضور 
داشته باشد. مقدار C‏ در واقع کنتراست بین اوزان است. مقدار صفر Of‏ بیان گر عدم وجود رابطه بین فاکتور مستقل با 
فاکتور وابسته (فرونشست) است. مقدار مثبت و بزرگتر نشان‌دهنده همبستگی مثبت و بیشتر بین فاکتور مربوطه با 
ظهور یا وقوع فرونشست است. در روش وزن‌دهی شواهد. بعد از محاسبه مقدار WE‏ در تمام طبقات فاکتورهاء ابتدا 
مقدار مجموع وزن‌های مثبت برای هر سلول روی نقشه به‌دست آمده و سپس احتمال نهایی محاسبه شد که مقدار آن 
بین صفر تا ۱ است. درنهایت با استفاده از نرم‌افزار ArcGIS‏ مدل به کار گرفته و نقشه میزان حساسیت فرونشست 
استخراج شد (عرفانیان و همکاران, ۱۳۹۵). 


] Pradhan & Lee 
2 Tweed 


Lal > ^Y‏ و مخاطرات محیطی شمارة چهارم 


۳- ارزیابی کارایی مدل‌های مورد استفاده 


ارزیابی دقت مدل‌ها در دو مرحله آموزش و اعتبارسنجی بر اساس روش منحنی مشخصه عامل گیرنده (ROC)‏ 
انجام شد. منحنی ROC‏ ب‌صورت نموداری ارائه می‌شود که در آن احتمال تشخیص صحیح نقاط حضور مدل با 
احتمال تشخیص صحیح نقاط عدم حضور مقایسه می‌شود (فیلپس " و همکاران ۲۰۰۹؛ الیت " و همکاران. ۲۰۱۱). 
سطح زیر نمودار (AUC)‏ به‌دست آمده به‌عنوان معیاری از قدرت تفکیک مدل در تشخیص نقاط حضور از عدم 
حضور موردتوجه قرار می‌گیرد. سطح زیر نمودار (AUC)‏ مدلی با قدرت پیشگویی بسیار کم» ۰/۵ است و مدلی 
کامل با قدرت پیشگویی بسیار ob;‏ دارای AUC‏ به میزان یک خواهد بود (فرامرزی و همکاران ۱۳۹۹). 


-٤‏ نتایج 
۱-۶6 عوامل مؤثر بر وقوع فرونشست 
نقشه عوامل مؤثر بر وقوع فرونشست زمین آبخوان دشت کاشان در شکل ۳ نشان داده شده‌اند. نقشه عامل فاصله 
از جاده نشان می‌دهد که تغییرات فاصله از جاده بین ۲۰۰ تا بیشتر از ۸۰۰ متر است و جاده‌ها در قسمت‌های شمال 
غربی» غرب و جنوب غربی آبخوان پراکنش بیشتری دارند (شکل (WILT‏ فاصله از معدن نیز بین ٤‏ تا ۲۰ کیلومتر 
متغیر است. معادن نیز بیشتر در قسمت‌های شمال غربی» غرب و جنوب غربی آبخوان و نزدیک به جاده‌ها پراکنش 
دارند (شکل ۳سب). میزان شیب در سطح دشت کمتر از ۲۱ درصد است و پراکندگی طبقات شیب در کل سطح 
آبخوان مشهود است. همچنین کمی بیش از AV‏ درصد آبخوان دشت LEIS‏ شیب کمتر از 1 درصد دارد (شکل Y‏ 
ج). میزان دامنه تغییرات انحنای طولی بین -۰/٩‏ الی ۰/٩‏ انحنای عرضی بین VA AVE‏ و انحنای سطح نیز بین 
۸- الی ۰/1۶ است و تغییرات نسبی این شاخص‌ها در کل سطح پراکنده است. با توجه به نتایج به‌دست آمده 
بیشتر سطح دشت در دامنه میانی قرار گرفته و درنتیجه, مسطح و بدون انحنا است (شکل‌های ۳د ۳سی CY‏ دامنه 
تغییرات ارتفاع در منطقه مطالعاتی بین ۸۰۳ الی ۱۱۷۱ متر است و از شمال و شمال غربی منطقه مطالعاتی به سمت 
جنوب و جنوب غربی ارتفاع از سطح دریا در حال افزايش است. همچنین حدود ۷۰ درصد از وسعت آبخوان دشت 
کاشان در دامنه ارتفاعی ۸۰۰ تا ۱۰۰0 متری قرار دارد (شکل PEAT‏ همچنین دامنه تغییرات سطح ایستابی در یک 
دوره ۱۲ ساله بین ۱۲/۰۸ YOUVE‏ متر است و تغییرات مذکور در جنوب غرب منطقه مطالعاتی بیشتر از سایر 
مناطق است. بیشترین تغییرات یعنی YOUN‏ متر نیز در منطقه حزاق است (شکل ۳ست). فاصله از گسل‌های منطقه 
از ۵۰۰ تا ۲۰۰۰ ga‏ متغیر است و پراکندگی گسل‌ها در شمال غربی. غرب. جنوب غربی و قسمتی از جنوب دشت 
بیشتر است (EY JE)‏ تنوع کاربری اراضی و سنگ‌شناسی نیز در منطقه مطالعاتی بالا است. به‌طوری که منطقه 


] Phillips 
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سال يازدهم ارزیابی کارایی مدل‌های دو متغیره در Oe‏ حساسیت پذیری AY se‏ 


مطالعاتی دارای ۸ نوع کاربری و ٩‏ نوع لیتولوژی است.کاربری‌های مرتع» کشاورزی و اراضی GLP cie‏ بیشترین 
سطح کاربری‌ها هستند و در میان آن‌ها Ge‏ کاربری مرتع بیشترین وسعت منطقه را به ميزان ۲ درصد به خود 
احتصاص داده است. واحد QU.‏ (تراس‌های قدیمی) در بخش لیتولوژی نیز دارای بیشترین وسعت (۵۳/۳۹ درصد) 
در منطقه مطالعاتی است. pls‏ یافته‌ها نیز نشان داد که حدود M‏ درصد منطقه مطالعاتی متشکل از تراس‌های قدیمی, 
ورق شن و ماسه و فلات رسی است (شکل‌های ۳-ح و ۳-س). تراکم چاه‌های بهره‌برداری نیز بین ۰/۳ تا ۷ تعداد 
در کیلومترمربع متغیر است و پراکندگی چاه‌ها در مرکز و جنوب شرقی دشت بیشتر از سایر مناطق است (شکل -Y‏ 
ش). فاصله از آبراهه‌ها بین ۲۰۰ تا ۸۰۰ متر متغیر است و الگوی شبکه آبراهه‌ها در منطقه مطالعاتی به‌صورت شاخه 


درختی است و دامنه تغییرات پایین OUS‏ از تعدد آبراهه‌ها در سطح دشت است (شکل LY‏ 


شمارة چهارم 
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- چیا مرس که 
شکل ۳- عوامل مورداستفاده در تعیین مناطق دارای فرونشست 


-۲-٤‏ وزن‌دهی طبقات عوامل مؤثر با استفاده از مدل‌های دو متغیره 
حدول A‏ وزن‌های ae‏ محاسبه شده cub‏ عوامل مؤثر بر فرونشست زمین بر اساس مدل‌های مورداستفاده در 
منطقه مطالعاتی را نشان می‌دهد. در این جدول همچنین مساحت هر طبقه برحسب هکتار و تعداد نقاط دارای 
فرونشست به تفکیک هر عامل نشان داده شده است. بر اساس نتایج به‌دست آمده بیشترین تعداد نقاط فرونشست 


M‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


مشاهداتی مورداستفاده در مرحله آموزش, در جهات جغرافیایی جنوب و شمال غربی» فاصله بیشت از ۰ متری 
گسل‌هاء شیب‌های کمتر از ۱۰/۵ درصد و همچنین ارتفاع کمتر از ۱۰۰۰ متری از سطح دریا واقع شده‌اند. همچنین 
«JS‏ فروند 


ت‌های ثبت شده در کاربری‌های اراضی nb‏ اراضی کشاورزی, مرتع و تنها یک مورد در منطقه 
مسکونی مشاهده شده است. بر اساس نتایج انحنای عرضی» انحنای طولی و انحنای ترکیبی» بیشترین میزان پراکنش 
فرونشست‌ها در مناطق مقعر دشت مشاهده شده است. تغییرات سطح ایستابی نیز نشان می‌دهد که در کلیه طبقات 
تغییر سطح ایستابی» فرونشست مشاهده شده است. بر اساس نتایج به‌دست آمده فرونشست‌ها بیشتر در مناطق دارای 
تراکم چاه کمتر رخ داده است. همچنین IS‏ فرونشست‌های ثبت شده در واحدهای سنگ‌شناسی Qe QCL‏ و Qt‏ 
ool NEAL laica‏ ۲۰۰-۰ و 1۱۰-۲۰۰ متری از آبراهه‌ها و نیز در 
فاصله بیشتر از ۸۰۰ متری از جاده‌ها بیشترین فرونشست اتفاق افتاده است. همچنین در کلیه طبقات فاصله از معادن 


فرونشست co‏ داده است. 


جدول ۱- وزن‌های نهایی روش‌های دمپستر- شفر. وزن‌دهی شواهد. انترویی شانون و نسبت فراوانی 
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سال یازدهم ارزیابی کارایی مدل‌های دو متغیره در تعیین حساسیت پذیری AV e‏ 
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فاصله از آبراهه 


فاصله از جاده 


سال يازدهم ارزیابی کارایی مدل‌های دو متغیره در Oe‏ حساسیت پذیری ^A Kis‏ 


۳-۶- نقشه‌های حساسیت فرونشست حاصل از مدل‌های مختلف 

نقشه حساسیت‌پذیری فرونشست منطقه مطالعاتی بر اساس مدل‌های فوق‌الذکر در شکل OUS ٤‏ داده شده است. 
در روش آنتروپی شانون» کلاس با حساسیت‌پذیری زیاد و خیلی 265 در مناطق شمال و شمال شرق دشت متمرکز 
هستند (شکل ٤-الف).‏ بر اساس نقشه حاصل از روش وزن شواهد. cle US‏ حساسیت‌پذیری زیاد و خیلی 265 
بیشتر در شمال غربی. غرب و جنوب شرقی دشت متمرکز هستند (شکل ۶-ب). همچنین نقشه حساسیت‌پذیری 
روش دمیستر- شفر نشان داد که بخش بسیار کمی از وسعت منطقه مطالعاتی دارای کلاس حساسیت‌پذیری زیاد و 
خیلی زیاد در شمال غربی. غرب و بخشی از جنوب و جنوب شرق دشت است CER)‏ بر اساس روش نسبت 
فراوانی» cla Us‏ با حساسیت‌پذیری زیاد و خیلی زیاد در شمال شرقی و شرق منطقه مطالعاتی متمرکز هستند 
(شکل (Gc‏ 

جدول Y‏ نیز مساحت (هکتار) و درصد هر یک از طبقات نقشه حساسیت‌پذیری فرونشست (بر اساس مرز نهایی 
در نظر گرفته شده) را به تفکیک چهار روش مورداستفاده نشان می‌دهد. بر اساس نتایج به‌دست آمده طبقه با حطر 
خیلی زیاد در همه روش‌ها کمترین میزان مساحت را به خود اختصاص داده است و دوطبقه خطر کم و متوسط 


جدول ۲- مساحت و درصد هر کدام از طبقات نقشه‌های حساسیت پذیری فرونشست با استفاده از مدل‌های مورد 


de ۲‏ کم کم متوسط زیاد خیلی زیاد 
روس 
مساحت درصد مساحت درصد مساحت درصد مساحت درصد مساحت درصد 
وزن شواهد VAL YoYAA/Y YVo£t £WOWAN ¥Y/44 VSSEAN YV/AA 11 10/71 1/0 YODAT‏ 


Y/^ Wet ۱۳/۳ To yay UVY YY/A VY4/0 YVI £VAWN آنتروپی شانون‎ 


YA OVYV/A EN 47/۱ JA] ¿NUM 0۰/۲ (۱۱۹ WA oYoY/o دمپسترسشفر‎ 


۳/۹ “OWN ۱۳/۳ 4/4 Yv/o MEZ Ya (0-۵ WN "S ÉAY/Y فراوانی نسبی‎ 


جغرافیا و مخاطرات محیطی 
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-٤ Js‏ نقشه‌های حساسیت پذیری فرونشست زمین در سطح دشت کاشان 
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-t-é‏ ارزیابی مدل‌ها 

در این مطالعه» ارزیابی نقشه‌های حاصل از مدل‌های مختلف با استفاده از روش منحنی ROC‏ برای مجموعه 
داده‌های مشاهداتی در دو مرحله مدل‌سازی و اعتبارسنجی انجام شد. شکل‌های ۵ و 1 منحنی‌های ROC‏ و مقادیر 
شاخحص AUC‏ را برای مدل‌های مورداستفاده در بخش آموزش و اعتبارسنجی مدل نشان می‌دهد. بر اساس شکل ۵ 
در مرحله آموزش» مدل آنترویی شانون با میزان AUC‏ برابر با ۰/۸7۶ بهترین عملکرد و مدل وزن‌دهی شواهد با 
مقدار AUC‏ برابر با ۰/۱۳۰ پایین‌ترین کارایی را نسبت به سایر روش‌ها داشته‌اند. دو مدل نسبت فراوانی و تابع باور 
شهودی (دمپستر- شفر) به‌ترتیب با AUC polis‏ برابر با ۰/۸۶۰ و ۰/۷۷۵ در اولویت‌های دوم و سوم قرار داشتند. 
نتایج حاصل از ارزیابی نقشه‌های حاصل از مدل‌ها بر اساس داده‌های اعتبارسنجی در شکل ٦‏ نیز حاکی از آن بود که 
روش نسبت فراوانی با مقدار AUC‏ برابر با ۰/۸۹۸ دقت بالاتری نسبت به سایر روش‌ها داشته است. پس‌ازآن, 
مدل‌های آنتروپی شانون (۰/۸1۶- (AUC‏ تابع باور شهودی (۰/۷۸۹- (AUC‏ و وزن‌دهی شواهد (AUC -*/WA)‏ 


در رده‌های بعدی از نظر کارایی قرار گرفتند. 


_— 
| J —— Hi if 

g 3 
H t 

€ 
3 4 int 1 3 a 

BANC EM 
na 7 A ۴9 
1 - Specificity 1 - Specificity 
سنحی مدل‌های‎ uel در بخش‎ AUC مقادیر‎ A در بخش آموزش مدل‌های شکل‎ AUC polis —0 شکل‎ 
مورداستفاده مورداستفاده‎ 


۵- بحث و نتیجه گیری 


فرونشسست می تواند در اثر پدیده‌های طبیعی مانند انحلال, تراکم نهشته‌ها و حرکات آرام پوسته و یا فعالیت‌های 


انسانی نظیر معدنکاری و برداشت آب‌های زیرزمینی ایجاد شود. در دمه‌های اخیر» فرونشست زمین همراه با صنعتی 
شدن» شهرنشینی و پمپاژ بی‌رویه از آب‌های زیرزمینی برای مصارف کشاورزی مشاهده شده است. در هنگام افت 
سطح آب زیرزمینی؛ زهکشی آب سبب کاهش فشار هیدرواستاتیکی در AU‏ رس و تحکیم و فشردگی این لایه‌ها و 
نشست زمین می شود. نتایج این پژوهش نشان داد که بیشترین میزان فرونشست‌ها در جهات جغرافیایی جنوب و 
شمال غربی واقع شده‌اند. در ple‏ جهات جغرافیایی نیز فرونشست در سطح منطقه مطالعاتی مشاهده شده است که 
اهمیت پایین جهت را در حساسیت فرونشست نشان می‌دهد. مختاری و همکاران (۱۳۹۸) و محمدی و همکاران 


(۲۰۱۹) نیز به اهمیت پایین جهت جغرافیایی و تأثیر Ol‏ بر فرونث 5 اشاره کرده‌اند. بر اساس نقشه فاصله از 


da LS‏ بیشترین میزان فرونشست‌ها در فاصله بیش از ۲۰۰۰ متری واقع شده‌اند که نشان‌دهنده تأثیر پایین این لایه بر 


اهمیت پایین گسل بر روی فرونشست اشاره کرده‌اند. این در Se‏ است که مختاری و همکاران (۱۳۹۸) در تسوج» 
محمدی و همکاران در دشت سمنان (۲۰۱۹) و تین‌بویی و همکاران (۲۰۱۸) در کره جنوبی به اثرات تشدید کننده 
گسل روی فرونشست اذعان کرده‌اند. بر اساس نقشه شیب دشت کاشان» تعداد 1۷ فرونشست در شیب کمتر از ۵/۷ 


درصد واقع شده‌اند. این موضوع نشان‌دهنده این است که در شیب‌های پایین‌تر خطر وقوع فرونشست بیشتر است 


Lal > ۹۲‏ و مخاطرات محیطی شمارة چهارم 


که مختاری و همکاران (YYAA)‏ و محمدی و همکاران (۲۰۱۹) نیز به این مورد اشاره کرده‌اند. دلیل این امر می‌تواند 
این باشد که در شیب‌های ol‏ امکان تشکیل آبخوان بیشتر است و لذا بر اثر برداشت آب در این مناطق» احتمال 


وقوع فرو 
اهمیت پایین‌تر شیب نسبت به سایر شاخص‌ها اشاره کرده‌اند. نتایج نشان داد که تعداد VV‏ مورد فرونشست در ارتفاع 


بین ۸۰۰ الی ۱۰۰۰ متر واقع شده‌اند و در ارتفاعات بالات ميزان وقوع فرو صفر بوده است که در تحقیقات 


محمدی و همکاران )14 (Y‏ اهمیت این شاخحص و اثرگذاری آن بر فرونشست بیان شده است. لایه کاربری اراضی 


نشان داد که تعداد ۵ مورد از à‏ 9 23 ت‌ها در اراضی pb‏ و کشاورزی ably‏ شده‌اند و این مورد نشان‌دهنده 


(۲۰۱۸) و محمدی و همکاران (۲۰۱۹) به این مهم اشاره کرده‌اند. یکی از علل این Ela y ey‏ منابع آب 
زیرزمینی در مناطق کشاورزی برای زراعت است که در پژوهش قرهچلو و همکاران (۰۰ع۱) به آن ASE‏ شده است. 
نتایج نشان داد که وقوع فرونشست در مناطقی رخ داده است که انحنای عرضی» طولی و انحنای ترکیبی دارای مقدار 


اهمیت پایین این شاخص‌ها اشاره کرده‌اند. بر اساس نقشه میزان تغییرات سطح ایستابی» تعداد فرونشست‌ها در 
مناطقی که میزان این تغییرات بین ۱ الی ۱۰ متر (کاهشی) بوده؛ بیشسترین مقدار بوده است. افت ۱۰ متری آب 
زیرزمینی همراه با تراکم چاه‌ها در US‏ شمال شرق و جنوب شرق دشت باعث فرونشست بیشتر در این مناطق شده 
است. تغییرات سطح آب زیرزمینی یک از مهم‌ترین شاخص‌های تأثیرگذار بر فرونشست منطقه مطالعاتی است که 
اهمیت این لایه در تحقیقات مختاری و همکاران (۱۳۹۸» جانباز فوتمی و همکاران (YS)‏ قربان‌زاده و همکاران 
(CY)‏ محمدی و همکاران (Y)‏ تین‌بویی و همکاران (Y)‏ رحمتی و همکاران (۲۰۱۹2)» کومار" و 
همکاران (۲۰۲۲) و نگویان" و همکاران (۲۰۲۲) نیز تأیید شده است. بر اساس نقشه تراکم چاه‌هاء تعداد وقوع 
فرونشست بیشتر در مناطقی رخ داده است که تراکم چاه کمتر از ۱/۳ کیلومتر در کیلومترمربع بوده است که این نتیجه 
بر خلاف انتظار بوده است. یکی از دلایل این امر می‌تواند برداشت آب بیشتر از چاه‌هایی که در مناطق با تراکم کمتر 
چاه قرار گرفته‌انده باشد. فرونشست‌های واقع در eU‏ لیتولوژی نشان داد که کل فرونشست‌های موجود در سطح 
آبخوان دشت در مناطق فلات رسی» ورقه‌های شن و ماسه و تراس‌های قدیمی مشاهده می‌شسوند. این نتایج 
نشان‌دهنده اثرات واحدهای مختلف لیتولوژی و اهمیت این شاخص در بحث فرونشست زمین است که در راستای 
تحقیقات رحمتی و همکاران (۲۰۱۹۵) و تین‌بویی و همکاران (۲۰۱۸) به‌عنوان یک شاخص مهم و تأثیرگذار است. 
پر اساس a‏ فاصله از آبراهه‌هاء تعداد 7۵ فرونشست در فاصله ۰ تا 1۰۰ متر مشاهده شد. در فواصل مختلف آبراهه 


] Kumar 
2 Nguyen 


سال eh‏ ارزیابی کارایی مدل‌های دو متغیره در Oe‏ حساسیت پذیری Meus‏ ۹۳ 


فرونشست‌هایی مشاهده شده است که OLE‏ از اهمیت پایین این شاخص و تأثیرگذاری کم آن روی فرونشست است 
که با تحقیقات مختاری و همکاران (۱۳۹۸) و رحمتی و همکاران (Y*34)‏ در یک راستاست. a‏ فاصله از حاده‌ها 
نشان داد که بیشترین میزان وقوع فرونشست در فاصله بیشتر از ۰ متر رخ داده است و نشان‌دهنده آهمیت oak‏ 
این pa US‏ بر روی فرونشست زمین است که در تحفیقات محمدی و همکاران (۲۰۱۹) به آن تا کید شده است. تأیه 
فاصله از معادن نشان داد که بیشترین ميزان وقوع فرونشست در فاصله بین TIL‏ ۰ متر رخ داده است. لازم 
به ذکر است که در همه فاصه‌ها از معادن موجود. فرونشست‌هایی مشاهده شده است که نشان‌دهنده اثر پایین این 

همچنین نتایج حاصل از ارزیابی نقشه‌های حاصل از مدل‌های دومتغیره بر اساس داده‌های اعتبارسنجی در منطقه 
است. «oU «s‏ مدل‌های آنترویی شانون» تابع باور شهودی و وزن‌دهی شواهد در رده‌های Sir‏ از نظر کارایی قرار 
فراوانی نتایج بهتری را ارائه کرده است. 

سایر یافته‌ها نشان داد که polis‏ کارایی مدل‌های دومتغیره در مرحله آموزش مدل بیشتر از مرحله اعتبارسنجی 
اعتبارسنجی تقسیم می‌شوند. در مرحله آموزش, نقشه حساسیت‌پذیری تهیه شده بر اساس حضور bu‏ فرونشست 
در طبقات عوامل به‌دست می‌آید. اما نقاط فرونشست مشاهداتی در مرحله اعتبارسنجی بدون هیچ دخل و تصرفی با 
نقشه نهایی حساسیت‌پذیری قطع داده می‌شوند. لذا این yl‏ قابل‌انتظار است که دقت در مرحله آموزش بالاتر از مرحله 
اعتبارسنجی باشد. 
اجرای برنامه‌های مختلف برای کاهش خطرات مرتبط با وقوع daz‏ فرونشست زمین استفاده شود. 
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